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t^ ' Abstract 
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\^ . We consider the Maldacena conjecture applied to the near horizon geometry of a Dl- 



brane in the supergravity approximation and consider the possibility of testing the 

0^ ■ conjecture against the boundary field theory calculation using DLCQ. We propose the 

two point function of the stress energy tensor as a convenient quantity that may be 

Q^' computed on both sides of the correspondence. On the supergravity side, we may invoke 

^ ■ the methods of Gubser, Klebanov, Polyakov, and Witten. On the field theory side, we 

derive an explicit expression for the two point function in terms of data that may be 

J^ ' extracted from a DLCQ calculation at a given harmonic resolution. This gives rise to a 

;_^ ■ well defined numerical algorithm for computing the two point function, which we test in 

the context of free fermions and the 't Hooft model. For the supersymmetric Yang-Mills 

theory with 16 supercharges that arises in the Maldacena conjecture, the algorithm 

is perfectly well defined, although the size of the numerical computation grows too 

fast to admit any detailed analysis at present, and our results are only preliminary. 

We are, however, able to present more detailed results on the supersymmetric DLCQ 

computation of the stress energy tensor correlators for two dimensional Yang Mills 

theories with (1,1) and (2,2) supersymmetries. 
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1 Introduction 

There has been a great deal of excitement during this past year following the realization 
that certain field theories admit concrete realizations as a string theory on a particular 
background [1]. By now many examples of this type of correspondence for field theories 
in various dimensions with various field contents have been reported in the literature (for 
a comprehensive review and list of references, see [2]). However, attempts to apply these 
correspondences to study the details of these theories have only met with limited success so 
far. The problem stems from the fact that our understanding of both sides of the correspon- 
dence is limited. On the field theory side, most of what we know comes from perturbation 
theory where we assume that the coupling is weak. On the string theory side, most of 
what we know comes from the supergravity approximation where the curvature is small. 
There are no known situations where both approximations are simultaneously valid. At the 
present time, comparisons between the dual gauge/string theories have been restricted to 
either qualitative issues or quantities constrained by symmetry. Any improvement in our 
understanding of field theories beyond perturbation theory or string theories beyond the 
supergravity approximation is therefore a welcome development. 

In this note we raise the Supersymmetric Discrete Light Cone Quantization (SDLCQ) of 
field theories [3, 4, 5] to the challenge of providing quantitative data which can be compared 
against the supergravity approximation on the string theory side of the correspondence. We 
will work in two space-time dimensions where the SDLCQ approach provides a natural non- 
perturbative solution to the theory. In general, attempts to improve the field theory side 
beyond perturbation theory seem like a promising approach in two space-time dimensions 
where a great deal is already known about field theories beyond perturbation theory. 

We will study the field theory/string theory correspondence motivated by considering 
the near-horizon decoupling limit of a Dl-brane in type IIB string theory [6]. The gauge 
theory corresponding to this theory is the Yang-Mills theory in two dimensions with 16 
supercharges. Its SDLCQ formulation was recently reported in [7]. This is probably the 
simplest known example of a field theory/string theory correspondence involving a field 
theory in two dimensions with a concrete Lagrangian formulation. 

A convenient quantity that can be computed on both sides of the correspondence is the 
correlation function of gauge invariant operators [8, 9]. We will focus on two point functions 
of the stress-energy tensor. This turns out to be a very convenient quantity to compute for 
many reasons that we will explain along the way. Some aspects of this as it pertains to a 
consideration of black hole entropy was recently discussed in [10]. There are other physical 
quantities often reported in the literature. In the DLCQ literature, the spectrum of hadrons 



is often reported. This would be fine for theories in a confining phase. However, we expect 
the SYM in two dimension to flow to a non-trivial conformal fixed point in the infra-red 
[6, 11]. The spectrum of states will therefore form a continuum and will be cumbersome to 
handle. On the string theory side, entropy density [12] and the quark anti-quark potential 
[12, 13, 14] are frequently reported. The definition of entropy density requires that we 
place the field theory in a space-like box which seems incommensurate with the discretized 
light cone. Similarly, a static quark anti-quark configuration does not fit very well inside a 
discretized light-cone geometry. The correlation function of point-like operators do not suffer 
from these problems. We should mention that there exists interesting work on computing 
the QCD string tension [15, 16] directly in the field theory. These authors find that the QCD 
string tension vanishes in the supersymmetric theories which is consistent with the power 
law quark anti-quark potential found on the supergravity side. 

2 Correlation functions from supergravity 

Let us begin by reviewing the computation of the correlation function of stress energy tensors 
on the string theory side using the supergravity approximation. The computation is essen- 
tially a generalization of [8, 9]. The main conclusion on the supergravity side was reported 
recently in [10] but we will elaborate further on the details. The near horizon geometry of a 
Dl-brane in string frame takes the form 
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In order to compute the two point function, we need to know the action for the diagonal 
fluctuations around this background to the quadratic order. What we need is an analogue 
of [17] for this background which unfortunately is not currently available in the literature. 
Fortunately, some diagonal fluctuating degrees of freedom can be identified by following the 
early work on black hole absorption cross-sections [18, 19]. In particular, we can show that 
the fluctuations parameterized according to 

ds"" = (l + fix'', U) + g{x\ U)) gooidxy + (l + 5/(xO, U) + g{x\ U)) gu{dx^f 
+ (l + f{x\ U) + g{x\ U)) guudU' + (l + /(x°, U) - ^^(x°, U)^ gnnd^'j 
e^ = (l + 3/(a;°,f/)-^(a;°,f/))e'^« (2) 



will satisfy the equations of motion 

9"(U)^^9'iU)-§,iU)-'-^%I^,m ^ (3) 

by direct substitution into the equations of motion in 10 dimensions. We have assumed 
without loss of generality that these fluctuation vary only along the x^ direction of the world 
volume coordinates like a plane wave e*^"^ . The fields f{U) and g{U) are scalars when the 
Dl-brane is viewed as a black hole in 9 dimensions; in fact there are the minimal and the 
fixed scalars in this black hole geometry. In 10 dimensions, however, we see that they are 
really part of the gravitational fluctuation. We expect therefore that they are associated 
with the stress-energy tensor in the operator field correspondence of [8, 9]. In the case of 
the correspondence between A/" = 4 SYM and AdS^ x S5, superconformal invariance allowed 
the identification of operators and fields in short multiplets [20]. For the Dl-brane, we do 
not have superconformal invariance and this technique is not applicable. In fact, we expect 
all fields of the theory consistent with the symmetry of a given operator to mix. The large 
distance behavior should then be dominated by the contribution with the longest range. The 
field f{k^, U) appears to be the one with the longest range since it is the lightest field. 

The equation (3) for f{U) can be solved explicitly in terms of the Bessel's function 



fiU) = U-'K,„{^lQTi^gl^NU-'k). (4) 

By thinking of f{U) in direct analogy with the minimally coupled scalar as was done in 
[8, 9], we can compute the fiux factor 
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up to a numerical coefficient of order one which we have suppressed. We see that the leading 
non-analytic (in k"^) contribution is due to the k^ term, whose Fourier transform scales 
according to^ 

{0{x)Om = ^^. (6) 

This result passes the following important consistency test. The SYM in 2 dimensions with 
16 supercharges have conformal fixed points in both UV and IR with central charges of 
order A^^ and A^, respectively. Therefore, we expect the two point function of stress energy 
tensors to scale like N'^/x'^ and N/x'^ in the deep UV and IR, respectively. According to the 

7-p - '^'^-P) 19 + 2p-p^ 

It is not difficult to show that for a generic p-brane, {O{x)O{0)) — N ^-v gy^ '' x ^-p 



analysis of [6] , we expect to deviate from these conformal behavior and cross over to a regime 
where supergravity calculation can be trusted. The cross over occurs at a; = 1/ gYMvN and 
X = vN/gYM- At these points, the A^ scaling of (6) and the conformal result match in the 
sense of the correspondence principle [21]. 

3 Correlation functions from DLCQ 

The challenge then is to attempt to reproduce the scaling relation (6), fix the numerical 
coefficient, and determine the detail of the cross-over behavior using SDLCQ. Ever since 
the original proposal [22], the question of equivalence between quantizing on a light-cone 
and on a space-like slice have been discussed extensively. This question is especially critical 
whenever a massless particle or a zero- mode in the quantization is present. It is generally 
believed that the massless theories can be described on the light-cone as long as we take 
m ^ as a limit. The issue of zero mode have been examined by many authors. Some 
recent accounts can be found in [23, 24, 25, 26, 27]. Generally speaking, supersymmetry 
seems to save SDLCQ from complicated zero-mode issues. We will not contribute much to 
these discussions. Instead, we will formulate the computation of the correlation function of 
stress energy tensor in naive DLCQ. To check that these results are sensible, we will first do 
the computation for the free fermions and the 't Hooft model. Extension to SYM with 16 
supercharges will be essentially straightforward, except for one caveat. In order to actually 
evaluate the correlation functions, we must resort to numerical analysis at the last stage of 
the computation. For the SYM with 16 supercharges, this problem grows too big too fast 
to be practical on desk top computer where the current calculations were performed. We 
can only provide an algorithm, which, when executed on an much more powerful computer, 
should reproduce (6). Nonetheless, the fact that we can define a concrete algorithm seems to 
be a progress in the right direction. One potential pit-fall is the fact that the computation 
may not show any sign of convergence. If this is the case, or if it converges to a result at 
odds with (6), we must go back and re-examine the issue of equivalence of forms and the 
issue of zero modes. 

The technique of DLCQ is reviewed by many authors [5, 28] so we will be brief here. The 
basic idea of light-cone quantization is to parameterize the space using light cone coordinates 
x'^ and x~ and to quantize the theory making x~^ play the role of time. In the discrete light 
cone approach, we require the momentum p_ = p~^ along the x~ direction to take on discrete 
values in units of p~^ /k where p"*" is the conserved total momentum of the system and k is an 
integer commonly referred to as the harmonic resolution. One can think of this discretization 
as a consequence of compactifying the x^ coordinate on a circle with a period 2L = 27ik/p'^. 



The advantage of discretizing the hght cone is the fact that the dimension of the Hilbert space 
becomes finite. Therefore, the Hamihonian is a finite dimensional matrix and its dynamics 
can be solved explicitly. In SDLCQ one makes the DLCQ approximation to the supercharges 
and these discrete representations satisfy the supersymmetry algebra. Therefore SDLCQ 
enjoys the improved renormalization properties of supersymmetric theories. Of course, to 
recover the continuum result, we must send k to infinity and as luck would have it, we find 
that SDLCQ usually converges faster than the naive DLCQ. Of course, in the process, the 
size of the matrices will grow, making the computation harder and harder. 

Let us now return to the problem at hand. We would like to compute a general expression 
of the form 

F{x-,x+) = {O{x-,x+)O{0, 0)) . (7) 

In DLCQ, where we fix the total momentum in the x~ direction, it is more natural to 
compute its Fourier transform 

F{P-,x+) = ^(O(P_,x+)O(-P_,0)) . (8) 

This can naturally be expressed in a spectrally decomposed form 



3.1 Free Dirac Fermions 

Let us first consider evaluating this expression for the stress-energy tensor in the theory of 
free Dirac fermions as a simple example. The Lagrangian for this theory is 

£ = ii|r^ ^ _ ^i|r^ (10) 

where for concreteness, we take 7" = 0-^,7"^ = icr^ and we take \E' = 2~^^^{'^). In terms of 

A 

the spinor components, the Lagrangian takes the form 

T/TTl 

Since we treat x^ as time and since x does not have any derivatives with respect to x^ in the 

Lagrangian, it can be eliminated from the equation of motion, leaving a Lagrangian which 

depends only on ip: 

Tn 1 

c = trd+^+t—rQ-^- (12) 



We can therefore express the canonical momentum and energy as 

P = dx^ iilj*d-ilj 

/im 1 

dx- - — ^Y^ • (13) 

In DLCQ, we compactify x^ to have period 2L. We can then expand ip and ip* in modes 
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^* = -L (6(-n)e'T^"" + ci(n)e-'T^"") . (14) 



Operators b{n) and d{n) with positive and negative n are interpreted as a destruction and 
creation operators, respectively. In a theory with only fermions, it is customary to take 
anti-periodic boundary condition in order to avoid zero-mode issues. Therefore, n will take 
on half- integer values^. They satisfy the anticommutation relation 

{b{n), 6(-m)} = {d{n), d{-m)} = 6n,m ■ (15) 

Now we are ready to evaluate (9) in DLCQ. As a simple and convenient choice, we take 

0{-k) = \l dx- {iij*d-ij - i{d-ij*)ij) e-"^^' . (16) 

which is the Fourier transform of the local expression for P_ with the total derivative con- 
tribution adjusted to make this operator Hermitian. Therefore, this should be thought of as 
the T^^ component of the stress energy tensor. For reasons that will become clear as we go 
on, this turns out to be one of the simplest things to compute. When acted on the vacuum, 
this operator creates a state 

T++(-A;)|0) = ^ ( - - n J h{-k + n)rf(-n)|0) . (17) 

Since the fermions in this theory are free, the plane wave states 

|n) = 6(-A; + n)rf(-n)|0) (18) 

constitute an eigenstate. The spectrum can easily be determined by commuting these oper- 
ators: 

( k k \ 
M^ln) = 2P_P+k) =mM- + - In) (19) 

\n k — n ) 

which is simply the discretized version of the spectrum of a two body continuum. All that 
we have to do now is calculate eigenstates of the actual theory we are interested in and 
^In SDLCQ one must use periodic boundary condition for all the fields to preserve the supersymmetry. 



to assemble these pieces into (9), but we can do a little more to make the result more 
presentable. The point is that since (9) is expressed in mixed momentum/position space 
notation in Minkowski space, the answer is inherently a complex quantity that is cumbersome 
to display. For the computation of two point function, however, we can go to position space 
by Fourier transforming with respect to the L variable. After Fourier transforming, it is 
straight forward to Euclideanize and display the two point function as a purely real function 
without loosing any information. To see how this works, let us write (9) in the form 
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The quantity inside the absolute value sign is designed to be independent of L. Now, to 
recover the position space form of the correlation function, we inverse Fourier transform 
with respect to P_ = /ctt/L. 
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The integral over L can be done explicitly and gives 
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where K4^{x) is the 4-th modified Bessel's function. We can now continue to Euclidean space 
by taking r^ = 2x~^x~ to be real and considering the quantity 
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This is a fundamental result which we will refer to a number of times in this paper. It 
has explicit dependence on the harmonic resolution parameter /c, but all dependence on 
unphysical quantities such as the size of the circle in the x~ direction and the momentum 
along that direction have been canceled. For the free fermion model, (23) evaluates to 
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with M^ given by (19). The large k limit can be gotten by replacing n —* kx and ^ J2n ~^ 
/q dx. We recover the identical result using Feynman rules. For r <^ m~^, this behaves like 
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3.2 't Hooft Model 

Let us now turn to a slightly more interesting problem of computing the correlation function 
of T++ in 't Hooft's model of two dimensional QCD [29] in the large N limit. This theory 
has two characteristic scales, one determined by the mass of the quarks and the other by 
the strength of the gauge coupling (yfy^iV. To a large extent, this is a solvable model. The 
spectrum and the wave function of the hadrons are encoded in a one parameter integral 
equation that can be handled in many ways. A thorough analysis of this model including 
the discussion of asymptotic behavior of certain correlation functions can be found in [30]. 
This is clearly a very mature subject. 

Applying DLCQ to the 't Hooft model is tantamount to placing 't Hooft's integral equa- 
tion for the meson spectrum on a lattice. The lattice is in the light-cone momentum space, 
which is precisely what is expected when the light-cone is compactified on a circle. The 
DLCQ of the 't Hooft model was analyzed in detail by [31, 32]. 

Our goal here is to show that the computation of (9) is straight forward and that it 
generates sensible answers. In fact, nothing could be simpler. The 't Hooft model is nothing 
more than a gauged version of the free fermion model. The Lagrangian for this theory is 
simply 

C = ^F^ + if'0^ - mf'^ . (26) 

^9ym 

We choose the light cone gauge A~^ = which is customary. One then finds that the A~ 

component of the gauge field is non-dynamical and can be eliminated using the equation 

of motion, just like the x component of the spinor in the free fermion model. Expressing 

everything in terms of the only dynamical field in this theory which is ip, one finds the 

canonical energy and momentum operators to take the form 

P = dx~ iilj*d-ilj 

'd-^l-'^rU-'-k'j-^A (27) 



where J_ = ipip*. All that changed in comparison to the free fermion model is the addition of 
a current exchange term in the light-cone Hamiltonian. Therefore, all we have to do here is 
to perform the identical computation specified by (23), but using the modified Hamiltonian, 
and letting M^ and \n) be the spectrum and the wavefunction of the n-th meson state in the 
spectrum. Since in DLCQ we are always working with a finite dimensional representation of 
the Hamiltonian dynamics, a small change in the form of the Hamiltonian matrix causes no 
particular difficulty. Let us discuss the result of such a computation. We will consider the 
case when Oym^ ^ ^^ so that the effect of the gauge interaction is strong. The spectrum 
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Figure 1: (a) Mass spectrum in units of gy^jN/n as a function of 1/k and (b) Log-Log plot 
of the Correlation function -^(T++(r)T++(0)) as a function of r for 't Hooft model computed 

using DLCQ with k = 5, k = 10, k = 50, and k = 100 in the units where gyMy Nj-n = 10^ 
and m = 10~^. The dotted line is for the free fermions with the same mass. 

can be computed reliably for large k and is in agreement with the results reported in [29] 
(see figure l.a). Since the state (17) created by operator (16) is odd under parity n ^^ k — n, 
only parity odd states contribute in the spectral decomposition. For r ^ l/gyAiVN, we 
expect the correlation function to behave just like the free fermion. The lightest meson in 
the parity odd sector has a mass of order gynVN. Due to the presence of this mass-gap, 
for r ^ l/gyMvN, we expect to see an exponential damping of the correlation function. 

This is precisely the behavior we seem to be finding. In figure l.b, we illustrate the result 
of computing (23) by first constructing the mass matrix M^ symbolically, then evaluating 
the spectrum and the eigenf unctions numerically and assembling the pieces. We have chosen 
to set gYMyN/ir = 10^ and m = 10~^. We have tried harmonic resolutions k = 5, k = 10, 
k = 50, and k = 100. Remarkably, the computation at a low harmonic resolution k seems 
not to be so far off from the result found using larger values of k. The correlation function 
appears to have more or less converged by the time we reach k = 100. We have also included 
a plot for gyM^ = with the same mass for comparison. 

One of the reasons why the convergence is relatively rapid is the fact that the matrix 
element {0\T^^ {k)\n) in (23) is not sensitive to the structure of the 't Hooft wave function 
at the boundaries. To see this more clearly, recall that in the continuum limit, we expect 
the eigenfunctions \n) to behave as x^ near x = where (3 is determined by 
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are effectively exclianging integral expression like 




Jo 


= 1 


by a discretized sum 
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whose leading correction for a > 1 is dominated by the terms of order 1/n. For a less than 
one, however, the leading correction is controlled by terms of order (-)". In computing the 
form factor for T"*"^(A;), we were fortunate to have only encountered an integral whose end 
point behavior went as x^^^. Had we instead chosen to compute two point function of a 
scalar operators like i/jiIj or ip'j^ip, we would have considered states 

mL /I 1 



j dx- ijij{x')e''ir-- \Q) = J2^(--T ) b{-n)d{-k + n)\0) 

/da;^^7^^(x^)e-'i^"nO) = Y — (- + t^^ K-n)d{-k + n)\Q) (31) 
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which gives rise to a pole near x = and x = 1 in the continuum limit. Therefore, the 
matrix element in the continuum limit behaves near x = as x^ . To control the error in 
this case, we must take 

5 « ' (32) 

or equivalently 

n > IX) (33) 

which grows exponentially with respect to f3~^. We would have had to work much harder if 
we had sent m/ gyuyN to zero. Luckily, this is not the case for with the T^"*" operator. 

3.3 Supersymmetric Yang-Mills theory with 16 supercharges 

Finally, let us turn to the problem of computing the two point function of the T^+ operator 
for the SYM with 16 supercharges. Just as in the 't Hooft model, adopting light-cone 
coordinates and choosing the light-cone gauge will eliminate the gauge boson and half of 
the fermion degrees of freedom. The most significant change comes from the fact that 
the fields in this theory are in the adjoint rather than the fundamental representations 
and the theory is supersymmetric. This does not cause any fundamental problem in the 
DLCQ formulation of these theories. Indeed, the SDLCQ formulation of this [7] as well 
as many other related models with adjoint fields have been studied in the literature. The 
main difficulty comes from the fact that in supersymmetric theories low mass states such as 

10 



tT[b{—ni)b{—n2)b{—k + rii + ^2)]|0) with an arbitrary number of excited quanta, or "bits," 
appear in the spectrum. This means that for a given harmonic resolution k, the dimension 
of the Hilbert space grows hke exp(v^), which is roughly the number of ways to partition k 
into sums of integers. 

The fact that the size of the problem grows very fast is somewhat discouraging from a 
numerical perspective. Nevertheless, it is interesting to note that DLCQ provides a well 
defined algorithm for computing a physical quantity like the two point function of T++ that 
can be compared with the prediction from supergravity. In the following, we will show that 
this can be computed for the SYM theory by a straight forward application of (23), just as 
we saw in the case of the 't Hooft model. 



The authors of [7] have shown that the momentum operator P^ is given by 
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The local Hermitian form of this operator is given by 
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where X and u are the physical adjoint scalars and fermions respectively, following the 
notation of [7]. When discretized, these operators have the mode expansion 
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Therefore, (L/7r)(0|T++(— A;)|n) is independent of L and can be substituted directly into (23) 
to give an explicit expression for the two point function. 

We see immediately that (37) has the correct small r behavior, for in that limit, (37) 
asymptotes to (assuming rif, = Uf) 
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which is what we expect for the theory of rib free bosons and Uf free fermions in the large k 
limit. 
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Figure 2: (a) The spectrum as a function of 1/k and (b) the Log-Log plot of the correlation 

^'^ '' v.s. r in the units where Qym^I'^ = 1 for A; = 3 



function (T++(x)T++(0)) (f+ 
and /c = 4. 
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Computing this quantity beyond the small r asymptotics, however, represents a formidable 
technical challenge. The authors of [7] were able to construct the mass matrix explicitly and 
compute the spectrum for k = 2, k = 3, and k = 4. Even for these modest values of the 
harmonic resolution, the dimension of the Hilbert space was as big as 256, 1632, and 29056 
respectively (the symmetries of the theory can be used to reduce the size of the calculation 
somewhat). In figure 2, we display results that parallel those obtained for the 't Hooft model 
earlier (figure 1) with the currently available values of k, except for the fact that we display 
the correlation function multiplied by a factor of 47r^r"^/A^^(2nf, + nj), so that it now asymp- 
totes to 1 (or in the logarithmic scale) in the k -^ oo limit. In this way any deviation from 
the asymptotic behavior 1/r^ is made more transparent. Note that with the values of the 
harmonic resolution k obtained at present, the spectrum in figure 2. a is far from resembling 
a dense continuum near M = 0. Clearly, we must probe much higher values of k before we 
can sensibly compare our field theory results with the prediction from supergravity. 

3.4 Supersymmetric Yang-Mills theory with less than 16 super- 
charges 

The computation of the correlator for the stress energy tensor in the (8,8) model is limited by 
our inability to carry out the computation for large enough harmonic resolution. It is the (8,8) 
model which we are ultimately interested in solving in order to compare against the prediction 
of Maldacena's conjecture in the supergravity limit. Nevertheless, the computation of the 
correlation function can just as well be applied to models with less supersymmetry. We will 
conclude by reporting the results of such a computation. 
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First, let us consider the theory with supercharges (1,1). This theory is argued not to 
exhibit dynamical supersymmetry breaking in [33, 34]. We can also provide a physicist's 
proof that supersymmetry is not spontaneously broken for this theory by adopting the argu- 
ment of Witten for the 2+1 dimensional SYM with Chern-Simons interaction [35]. In [35], 
the index of 2+1 dimensional SYM with gauge group SU{N) and 2 supercharges on i? x T^ 
was computed and was found to be non- vanishing for Chern-Simons coupling k^ > N/2. If 
the period L of one of the circles in T^ is sufficiently small, this theory is approximately 
the 2-dimensional SYM with (1,1) supersymmetry with gauge coupling (^1 = g'^/L and BF 
coupling k2 = k^L [36]. Imagine approaching this theory by taking L ^ keeping g2 and 
^3 fixed. In this limit, A;2 ^ in the units of g2 so the limiting theory is that of pure SYM 
with (1,1) supersymmetry and a vanishing BF coupling. Choosing different values of k^ cor- 
responds to a different choice in the path of approach to this limit. If we chose k^ > N/2, we 
are guaranteed to have a non-zero index for finite L. This means that there will be a state 
with zero mass in the L — > limit also, indicating that supersymmetry is not spontaneously 
broken in this limit. On the other hand, the index is not a well defined quantity in the L ^ 
limit, as a different choice of k^ will lead to a different value of the index in the L ^ limit. 
In fact, the index can be made arbitrarily large by taking k^ to be also arbitrarily large. 
This suggests that there are infinitely many states forming a continuum near m = 0. The 
index is therefore an ill defined quantity, akin to counting the number of exactly zero energy 
states on a periodic box as one takes the volume to infinity. 

This theory is also believed not to be confining [15, 16] and is therefore expected to 
exhibit non-trivial infra-red dynamics. 

The SDLCQ of the 1+1 dimensional model with (1,1) supersymmetry was solved in 
[37, 38], and we apply these results directly in order to compute (23). For simplicity, we 
work to leading order in the large A^ expansion. The spectrum of this theory for various 
values of k, and the subsequent computation of (23) is illustrated in figure 3. a. 

The spectrum of this theory at finite k, illustrated in figure 3. a, consists oi2k — 2 exactly 
massless states^, accompanied by large numbers of massive states separated by a gap. The 
gap appears to be closing in the limit of large k however. We have tried extrapolating the 
mass of the lightest massive state as a function of 1/k by performing a least square fit to a line 
and a parabola, giving the extrapolated value of M'^Ti/gyM^ = 1-7 and M'^n/gyM^ = —0.6, 
suggesting indeed that at large k, the gap is closed. This is consistent with the expectation 
that the spectrum is that of a continuum starting at M = discussed earlier, although one 
must be careful when the order of large N and large k limits are exchanged. At finite A^, we 
expect the degeneracy oi 2k — 2 exactly massless states to be broken, giving rise to precisely 
■^i.e. A: — 1 massless bosons, and their superpartners. 
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a continuum of states starting at M = as expected. 

In the computation of the correlation function illustrated in figure 3.b, we find a curious 
feature that it asymptotes to the inverse power law c/r'^ for large r. This behavior comes 
about due to the coupling (0|T"''+|ri) with exactly massless states \n). The contribution to 
(23) from strictly massless states are given by 
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T ^ 
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(39) 

We have computed this quantity as a function of 1/k and extrapolated to 1/k — > by fitting 
a line and a parabola to the computed values for finite 1/k. The result of this extrapolation 
is illustrated in figure 4. The data currently available suggests that the non-zero contribution 
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from these massless states persists in the large k limit. 

Let us now turn to the model with (2,2) supersymmetry. The SDLCQ version of this 
model was solved in [39]. The result of this computation can be apphed to (23). The result 
is summarized in figure 5. This model appears to exhibit the onset of a gapless continuum 
of states more rapidly than the (1,1) model as the harmonic resolution k is increased. Just 
as we found in the (1,1) model, this theory contains exactly massless states in the spectrum. 
These massless states appear to couple to T++|0) only for k even, and the overlap appears 
to be decreasing as k is increased. We believe that this model is likely to exhibit a power 
law behavior c/r'^ for 7 > 4 for the T"^+ correlator for r ^ QyaiVN in the large N limit. 
Unfortunately, the existing numerical data do not permit the reliable computation of the 
exponent 7. 



4 Conclusion 

In this article, we have provided a prescription for computing the correlation function of the 
stress energy tensor T^^ in the SDLCQ formalism, which may be readily compared with 
predictions provided by a supergravity analysis following the conjecture of Maldacena. Such 
a comparison requires non-perturbative methods on the field theory side, and the SDLCQ 
approach appears at first sight to be particularly well suited to this task. Unfortunately, 
at the present time, high enough resolution calculation have not been made to evaluate ex- 
pression (23) accurately enough in the case of SYM with 16 supercharges to reproduce (6). 
Significant progress is expected however when these calculation are moved from desk top 
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computers to a supercomputer. Ultimately the main obstacle will be that the number of 
allowed states in the SDLCQ wavefunctions grows exponentially with the resolution. Never- 
theless, a concrete well defined algorithm is a great starting point for further investigations, 
and additional insight may be gained by studying models with less supersymmetry, as we 
have done here. 

Thus, the answer to the question posed in the title is a qualified "yes." Equation (23) may 
be compared against (6) using SDLCQ, and unless other computational approaches such as 
the lattice technique catches up, the SDLCQ approach remains the most viable option. 
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